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Abstract. The foundations of gyrokinetic theory are reviewed with an 
emphasis on the apphcations of Lagrangian and Hamiltonian methods 
used in the derivation of nonlinear gyrokinetic Vlasov-MaxweU equa- 
tions. These reduced dynamical equations describe the turbulent evolu- 
tion of low-frequency electromagnetic fluctuations in nonuniform mag- 
netized plasmas with arbitrary magnetic geometry. 



1 Introduction 

The development of gyrokinetic theory was initially motivated by the need to 
describe complex plasma dynamics over time scales that are long compared to the 
short gyration time scale of a charged particle about a nonuniform magnetic field 
line. Thus, gyrokinetic theory was built upon a generalization of guiding-center 
theory (see Refs. and [2]). In Ref. (S], for example, Taylor showed that, while 
the guiding-center magnetic-moment invariant can be destroyed by low-frequency, 
short perpendicular-wavelength electrostatic fluctuations, a new magnetic-moment 
invariant can be constructed as an asymptotic expansion in powers of the amplitude 
of the perturbation field. This early result indicated that gyrokinetic theory was to 
be built upon an additional transformation beyond the guiding-center phase-space 
coordinates, thereby constructing new gyrocenter phase-space coordinates, which 
describe gyroangle-averaged perturbed guiding-center dynamics. The purpose of 
the present paper is to review the foundations of gyrokinetic theory by presenting 
the Lagrangian and Hamiltonian methods used in the derivation of self-consistent, 
energy-conserving gyrokinetic Vlasov-Maxwell equations describing the nonlinear 
turbulent evolution of low-frequency, short perpendicular-wavelength electromag- 
netic fluctuations in nonuniform magnetized plasmas. 

Although a comprehensive review of the applications of gyrokinetic theory mer- 
its a separate paper, we briefly point out, here, that applications are generally di- 
vided into applications based on linearized gyrokinetic equations and applications 
based on nonlinear gyrokinetic equations. On the one hand, linear gyrokinetic 
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theory (see Refs. |4]-|10|. for example) is concerned with the stabihty of inhomoge- 
neous magnetized plasmas perturbed by low-frequency electromagnetic perturba- 
tions (which preserve the gyrocenter magnetic moment, or first adiabatic invariant). 
Nonlinear gyrokinetic theory (see Refs. jll| - |16j V on the other hand, focuses its at- 
tention on the anomalous transport associated with low-frequency electromagnetic 
fluctuations in inhomogeneous magnetized plasmas. Gyrokinetic particle simulation 
techniques (see Refs. |17j-|20)l now play a major role in the investigation of low- 
frequency plasma turbulence and its associated transport in magnetized plasmas 



The remaining material is organized as follows. In Sec. [3 we introduce the 
variational principles in eight-dimensional extended phase space that are necessary 
to present a consistent time-dependent Hamiltonian perturbation theory for single- 
particle dynamics. Here, and throughout the paper, the language of differential 
forms is strongly emphasized within the context of Lagrangian and Hamiltonian 
mechanics 22 . In Sec. 13 the phase-space Lagrangian Lie-transform perturba- 
tion method, which allows a simultaneous dynamical reduction of the Hamiltonian 
and Poisson-bracket structure for single-particle dynamics, is briefly presented. In 
Sec. 01 the phase-space Lagrangian Lie-perturbation method is used to derive gyro- 
center Hamiltonian dynamics by a time-dependent phase-space transformation from 
guiding-center coordinates to gyrocenter coordinates. Here, we choose a Hamilton- 
ian formulation that preserves the guiding-center Poisson-bracket structure and 
carries all field perturbations onto the gyrocenter Hamiltonian. 

In Sec. [SJ we show how a self-consistent set of nonlinear gyrokinetic Vlasov- 
Maxwell equations can be derived from a reduced variational principle. The reduced 
variational principle is itself shown to be derived from a variational principle for 
the exact Vlasov-Maxwell equations also presented in Sec.|31 We also show how the 
Noether method can be used to derive exact conservation laws for the exact and 
gyrokinetic Vlasov-Maxwell equations. Lastly, we summarize the work presented 
here and discuss extensions of the standard gyrokinetic formalism in Sec. El 



This Section presents a brief introduction to the extended phase-space La- 
grangian formulation of charged-particle dynamics in a time-dependent electro- 
magnetic field. Here, the electromagnetic field is represented by the potentials 
(<&, A), while the eight-dimensional extended phase-space noncanonical coordinates 
2, — (x, V, w, t) include the position x of a charged particle (mass m and charge e), 
its velocity v, and the canonically-conjugate time-energy (t, w) coordinates. 

2.1 Single-particle Lagrangian dynamics in extended phase space. 

The phase-space Lagrangian (or Poincare-Cartan 22 ) one- form for a charged par- 
ticle in eight-dimensional extended phase space is expressed as 



r = [- K ^ mv] -dy. ~ w dt - Udr = V a{Z) dZ"- - n{Z) dr, (1) 



where summation over repeated indices is implied (latin letters a,b,c, ... go from 1 
to 8 while greek letters /i, i^, ... go from to 3), r denotes the Hamiltonian orbit 
parameter, and the extended phase-space Hamiltonian is 



m 



2 Variational Principles for 
Extended Single-Particle Phase-Space Dynamics 




— w. 



(2) 
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Here, H(z,t) denotes the standard time-dependent Hamiltonian, with z — (x, v), 
and the physical single-particle motion takes place on the subspace Ti-iZ) — H{z, t) — 
w = of extended phase space. 

Next, we introduce the single-particle action integral 

where Ta are known as the symplectic components of the phase-space Lagrangian 
r. The Principle of Least Action 6S ^ for single-particle motion in extended 
phase space 



= ST = SZ" 



dZ^ dZ^ J dZ'^ 
yields the extended phase-space Euler-Lagrange equations 



(3) 



dz" ^ fm\ _ dFA dz^ _ m_ 

dT - \dz- dz" dT ~ dz-' ^ ' 



where uJat denotes a component of the 8x8 Lagrange matrix u) 

The components of the inverse of the Lagrange matrix J = UJ~^, known as the 
Poisson matrix, are the fundamental Poisson brackets 

= {2'', Z'}z = J''\Z), (5) 

so that the Euler-Lagrange equations Q become the extended Hamilton's equations 

^ = r" §^ = (^". (6) 

The fundamental Poisson brackets J°''(Z) satisfy the Liouville identities (which 
follow from the incompressibility of the extended Hamiltonian flow) 

J^aiJJ^')=0, (7) 



where J{Z,) = \J det uj defines the Jacobian associated with the choice of eight- 
dimensional phase-space coordinates, and the Jacobi property 

' Jz-^^'^^'^-'^ 

obtained from the Jacobi identity 

{/, {5, h}} + {g, {h, /}} + {h, {/, g}} = 0, (8) 

which is valid for three arbitrary functions /, 5, and h. Note that, by using the 
identity cj^c J'^^ — S^, the Jacobi property may be rewritten as d w = 0, i.e., 

daUJbc + dbUJca + OcLOab = 0, 

which is always satisfied since w = dF is an exact two- form (i.e., ujab = da^b ~ 
SfcFa). Hence, any bracket derived through the sequence F— ^ijj = dr^J = lo~^ 
automatically satisfies the Jacobi identity ©. 

Lastly, using the extended phase-space Lagrangian |Q , the explicit form of the 
extended noncanonical Poisson bracket, defined by Eq. (jSJ, is 

df fdg 



dw V dt 



e dA 


dg\ 


(df 


e dA 


df\ 


dg 


mc dt 


d^l 


\ dt 


mc dt 


dv) 


dw 



I dg df ^\ eB df dg 

m \ dv dv I m^c av dv 
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Here, the noncanonical aspect of the Poisson bracket Q is exhibited by the ap- 
pearance of space-time derivatives of the vector potential A(x, t). The Hamiltonian 
dynamics in extended phase space is expressed in terms of Eqs. ©-© a,s 

dx _ 1 dn _ 
dr m 9v 

dv 1 e fdA dn dn B\ e / V „ 

— = +— —— + — X- =-E+-xB 

dr m mc \ dt dw dw m J m \ c 

dw _ dn e dA dn _ fd^ v dA 
dr dt mc dt 9v \dt c dt 

dt _ dn _ ^ 
dr dw 

Note that the Hamihonian orbit parameter r can be identified with the time coor- 
dinate t and, as expected, the energy coordinate is an invariant when the electro- 
magnetic fields are time- independent. 

2.2 Single-particle perturbation theory in extended phase space. It 

turns out that the eight-dimensional extended phase space provides a natural setting 
for time-dependent Hamiltonian perturbation theory. A variational formulation of 
single-particle perturbation theory can be introduced through the new phase-space 
Lagrangian one-form |24j 

r' = TadZ" ~ ndr ~ S de, (10) 

where the symplectic components Ta and the Hamiltonian n now depend on the 
perturbation parameter e and the scalar field S is the generating function for an 
infinitesimal canonical transformation that smoothly deforms a particle's extended 
phase-space orbit from a reference orbit (at e = 0) to a perturbed orbit (for e 7^ 0). 
From the phase-space Lagrangian H1U|) , we construct the action pat/i-integral — 
V evaluated along a fixed path C in the (r, e)-parameter space. 
The modified Principle of Least Action for perturbed single-particle motion in 
extended phase space 



= sr' = / 



dn , ( ds dTa. , , 

^•^^dZ - —dr - + ^]de 



(11) 



now yields the extended perturbed Hamilton's equations 
dZ°- 

= {z\n}^, (12) 



dr 

{Z\,S}^~'^{Z\Z-]^, (13) 



where Eq. H12|l is identical to Eq. lO except that the extended Hamiltonian and 
Poisson bracket now depend on the perturbation parameter e, while Eq. (|13l) de- 
termines how particle orbits evolve under the perturbation e-flow. 

We note that the order of time evolution (r-flow) and perturbation evolution 
(e-flow) should be immaterial (i.e., we may evolve the system along its reference 
orbit and then perturb it its final state, or perturb the initial state and evolve the 
perturbed orbit to the same final state). The commutativity of the two Hamiltonian 
(r, e) flows, therefore, leads to the path independence of the action integral / F' in 
the two-dimensional (r, e) orbit-parameter space. Thus, considering two arbitrary 
paths C and C with identical end points on the (r, e)-parameter space and calculat- 
ing the action path-integrals S'(j = V and Sy — J-^ F', the path-independence 
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condition (5^ = S'q) leads (by Stokes' Theorem for differential one- forms |25| ) to 



the condition 







r' 



r' = 



r' = 



dD 



dV, 



D 



where D is the area enclosed by the closed path dD = C — C, and the two-form 
dV on the (r, e)-parameter space is 



dV 



de A dr 



dZ'' 



dZ^ 

"3 — ^ab —5 — 

de dr 



deAdr {5, n\z 



m 



on on dz'' 
"aT ~^ 'dz^ "dT 

dTa 



OS 
dZ^ 



dFA dZ^ 
de ) dr 



de 



where Eqs. H12|l - H13|) were used. Hence, the condition of path independence requires 
that dr' = 0, which yields the Hamiltonian perturbation equation 



{5, - ^ - §^ > 



(14) 



relating the generating scalar field S to the perturbation-parameter dependence of 
the extended Hamiltonian {dtH) and Poisson bracket {dfTa)- Using the extended 
single-particle Hamiltonian (O, the Hamiltonian perturbation equation (|14|l be- 
comes 



9<i> e 9A 



(15) 



whose formal solution is 



'd7 



de 



c de 



where {d/dT)~^ denotes an integration along a perturbed Hamiltonian orbit. Here, 
the perturbed evolution operator d/dr is expanded in powers of e, with the lower- 
order operator do/dr considered to be explicitly integrable. In practice, the gener- 
ating function S is also expanded in powers of e: S = Si + e S2 + ■ ■ ■ , so that the 
first-order term is expressed as 



Si 



e $1 — e Ai 



vo 



(16) 



where vq = dox/dr — {x, 7io}o denotes the particle's unperturbed velocity. In or- 
der to determine the higher-order terms iS„ (for n > 2), a more systematic approach, 
based on applications of the Lie-transform perturbation method, is required. 



3 Lie-transform Perturbation Theory 

The Hamiltonian perturbation equation 1)14(1 arises naturally within the con- 
text of the dynamical reduction of single-particle Hamilton's equations p2(l through 
the elimination of fast orbital time scales. The most efficient method for deriving 
reduced Hamilton's equations is based on Hamiltonian |26| and phase-space La- 
grangian '27' Lie-transform perturbation methods. 
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3.1 Near-identity phase-space transformations. The process by which a 
fast time scale is removed from Hamilton's equations Ti.}z involves a near- 

identity transformation on extended particle phase space |28j : 

%: Z ^ Z{Z;e) = %Z, with Z(Z;0) = Z, (17) 

where e <C 1 denotes a dimensionless ordering parameter. Here, the near-identity 
transformation is explicitly expressed in terms of generating vector fields {Qi, C/2, ■■■)'■ 

Z\z,e) = Z^ + eQl + {g^, + i c;^^ g ^ + ... , (ig) 

where the nth-order generating vector field Qn is chosen to remove the fast time 
scale at order e" from the perturbed Hamiltonian dynamics. The new extended 
phase-space coordinates include the pair of fast action- angle coordinates (J, 9) and 
the reduced phase-space coordinates Zji such that the fast action J = jQ-\-e Ji + - ■ ■ 
is an exact invariant of the reduced Hamiltonian dynamics and the Hamiltonian dy- 
namics of the reduced coordinates Z_r is independent of the fast angle 9. The small 
dimensionless parameter can, therefore, be defined as e = (thJI)"^ <C 1, where t/j 
is the characteristic time scale of the reduced dynamics and {2-k/U) denotes the 
fast orbital time scale associated with the fast angle 9. 

Next, using the transformation H17|l . we define the push-forward operator on 
scalar fields |28i induced by the near-identity transformation H17|l : 

J-'^ : T y = T^^T, (19) 

i.e., T^^ transforms a scalar field J- on the phase space with coordinates Z into a 
scalar field T on the phase space with coordinates Z: 

J(Z) = T-^T(Z) ^ T{T-^Z) = T{Z). 

Since the transformation H17|) is invertible, i.e., there exists an inverse near- identity 
transformation 

: Z -> Z(Z;e) = T'^Z, with Z(Z;0) = Z, (20) 

we also define the pull-back operator j28| : 

l,:y-.T = TeJ, (21) 

i.e., Te transforms a scalar field T on the phase space with coordinates Z into a 
scalar field J- on the phase space with coordinates Z: 

T{z) = j,y{z) = T{%z) = y(z). 

Using the fact that the total r-derivative of a scalar field J- is itself a scalar field, 
we obtain an expression for the transformed operator d^/ dr defined as 

where { , ^-^ denotes the new transformed Poisson bracket and 

n = J-^n (23) 

denotes the transformed Hamiltonian. Once again, the new extended phase-space 
coordinates are chosen so that d^J/dr = {J, H}^ = and deZn/dr — {Zr, Ti.}^ 
be independent of the fast angle 9. The dynamical reduction of single-particle 
Hamiltonian dynamics consists in the construction of a fast invariant J with its 
canonically-conjugate fast-angle 9 becoming an ignorable coordinate. 
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3.2 Lie derivatives and Lie transforms. In Lie-transform perturbation 
theory j28| . the push- forward operator 1)19(1 is defined as 

= ••• exp(-e2L2) cxp(-eLi) (24) 

in terms of the nth-order Lie derivative L„ generated by the nth-order vector field 
On ■ In Eq. the Lie derivative L„H is defined as the scalar field 

LnH = g^daU. (25) 

The transformation of the Poisson bracket by Lie-transform methods, on the other 
hand, is performed through the transformation of the extended phase-space La- 
grangian, expressed as 

r = T^ip + d5, (26) 

where S denotes a (gauge) scalar field used to simplify the transformed Hamiltonian 
(|23|l . i.e., it has no impact on the new Poisson-bracket structure ZU = dP = d(T~^P) 
since d^5 = (i.e., d'^f^S — d'^^S = 0). In Eq. if^ . the nth-order Lie derivative L„P 
of a one-form P = P„ dZ'^ is defined as a one-form 



L„T = gi iuab dZ" + d {gi To) , (27) 

where, at each order, the terms d(- • • ) can be absorbed in the gauge term d5„. 

3.2.1 Transformed extended Poisson-bracket structure. We now write the ex- 
tended phase-space Lagrangian P = Pq -|- ePi and the extended Hamiltonian Ti. = 
Tio + cTi-i in terms of an unperturbed (zeroth-order) part and a perturbation (first- 
order) part. The Lie-transform relations associated with Eq. ((26|l are expressed (up 
to second order in e) as 

Poa = Poa, (28) 
Lla = Pla - g'l i^Oba + (29) 

P2a = - ^2 ^06a ^ 2 ^^^''"^ ~^ ^'^ba) + daS2, (30) 

where each scalar field 5„ (n > 1) is determined by requiring that the nth-order 
Hamiltonian 7i„ be independent of the fast orbital time scale. The simplest form 
for the new Poisson bracket { , j^" is obtained by choosing P = Pq = Pq, so that 
the condition P„ = (n > 1) yield a solution for the generating vector field t/„ 
expressed in terms of the scalar fields (5i, • • • ,iS„). 

For the first-order generating vector field 5i, the condition Pi = yields the 
following expression in terms of the scalar field Si : 

gi = Z'^lo -f Pib Jo^^ (31) 

where Jg^ denotes a component of the zeroth-order Poisson matrix. Next, for the 
second-order generating vector field C/2, the condition P2 = yields the following 
expression in terms of the scalar field 52: 

C/2 = {^2, Z^lg - - t/i CJif,c (32) 

where uJi \,c denotes the component of the first-order perturbed Lagrange matrix. 
The near-identity extended phase-space transformation l|18f) is, thus, expressed (up 
to first order in e) as 

= Z'' -f e({5i, Z'^lo + Pifc J^^-) + ••• , (33) 
and its explicit expression requires a solution of the scalar fields {Si ,•••)■ 
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3.2.2 Transformed extended Hamiltonian. By substituting the generating vec- 
tor fields 1)31(1 and H32|l into the Lie-transform relations associated with Eq. ((23|l : 

'Ho = ^Oj 

Hi = Hi - Q^daHo, 

^2 = - Gi daUo ~ \ {Hi + Hi) , 

we obtain the first-order and second-order terms in the transformed extended 
Hamiltonian: 

Hi ^ Hi — Via {2°, TYolo ^ i*^!' ^olo = A^i - {^i, Tiolo 7 ("^4) 
where /Ci = TYi — Fia {2^°, 7io}o denotes the effective first-order Hamiltonian, and 

H.2 = - {S2, H.o}q - 2(^"{'^i' H-i+Hi}^ + Tia {Z", Hi+Hi}^ 

+ {Z^ 7^0)0 • (35) 

The solution for the new first-order Hamiltonian H34|) is expressed in terms of the 
fast-angle averaging operation (• • • ) as 

Hi = (/Ci> = (Hi) - (Fi, {Z^ Ho}), (36) 

where the Poisson bracket { , } is the zeroth-order Poisson bracket { , }o (unless 
otherwise noted) and Si can be chosen such that {Si) = 0. The first-order scalar 
field Si is, thus, the solution to the perturbation equation 

{Si, Ho} = K.I = K.I - (/Ci) Si = {do/dT)-'iCi, (37) 

where {dQ/dr)^^ denotes an integration along an unperturbed extended Hamilton- 
ian orbit; note the similarity of Eq. 1)37(1 with Eq. III6I) . To lowest order in the 
fast orbital time scale, the unperturbed integration {do/dT)^^ICi = il^^ j ICidO 
involves an indefinite fast-angle integration. 

The solution for the new second-order Hamiltonian 1(35(1 can be simplified as 
follows. First, we introduce the Poisson-bracket identity (valid for four arbitrary 
functions /, g, h, and k) 

{/, 9} {h, k}-{f, h} {g, k} - {/, {g{h, k})} + {{g{h, /}), k} - g{h, {/, fc}}, 
which follows from the Jacobi identity (jHJ, so that the last term in Eq. 1(35(1 
G'^ uiab {Z\ Tio} = ({5i, FiJ {2", Ho} - {Si, Z"} {Fi„ Ho}) 

-I- ^l„({z^ Fib} {z", n] - {z^ z^} {Fi^, Ho}) 

can be rearranged and inserted back into Eq. ((35(1 to yield 

n2 = - {('^2 - ^ Fl,{z^ Si}^ , Hoj - 1^1, (^(/Ci) + ^ici 

- Via (^{Z^ ICi} + i Fib I 2^ {Z^ Ho}|) . (38) 

Since ^2 can also be chosen such that (52) = 0, the fast-angle average of the right 
side of Eq. ((38(1 yields the new second-order Hamiltonian 

'H2 = I (^Tia {Z", {Z^ Ho}} Fib^ - i ^ [Si, iCi}' 

~ (Fi, ({z^ Hi} ~ {z^ Fi4 {z^ ho})) . (39) 
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In the next Section, we will show that the first term in Eq. H39() corresponds to the 
expected quadratic nonlinearity associated with a perturbed canonical representa- 
tion in which only the perturbation vector potential Ai appears in the expression 
for the reduced momentum coordinates. The second term, on the other hand, 
corresponds to the low-frequency limit of the standard quadratic ponderomotive 
Hamiltonian |26j . Lastly, the third set of terms will be shown to vanish because 
the Poisson brackets {x, $i(x)} and {x, Ai(x)} both vanish. 

4 Nonlinear Lov^r-frequency Gyrocenter Hamiltonian Dynamics 

In this Section, we apply the Lie-transform perturbation methods to the dy- 
namical reduction of the perturbed dynamics of charged particles (mass m and 
charge e) moving in a background time-independent magnetic field Bp = V X Ag 
in the presence of low-frequency electromagnetic fluctuations represented by the 
perturbation four-potential A'^ — ($1, Ai), whose amplitude is ordered with a di- 
mensionless small parameter ^ 1. 

The eight-dimensional extended phase-space dynamics is expressed in terms 
of the extended phase-space Lagrangian F = Fq + £5 Fi, where Fq = [(e/c) Aq + 
TOv] ■ dx — w dt and Fi = (e/c) Ai • c?x, and the extended phase-space Hamiltonian 
H = Ho + esHi, where Ho = i''n/2) |vp — w and Hi = e^i. The extended 
Poisson bracket { , is obtained from the extended phase-space Lagrangian F by 
standard means j28l . as described earlier. Note that, while electrostatic fluctuations 
perturb the Hamiltonian alone, full electromagnetic fluctuations perturb both the 
Hamiltonian and the Poisson bracket. 

The standard gyrokinetic analysis for magnetized plasmas perturbed by low- 
frequency electromagnetic fluctuations |15j proceeds by a sequence of two near- 
identity phase-space transformations: a time-independent guiding- center phase- 
space transformation and a time-dependent gyrocenter phase-space transformation. 

4.1 Unperturbed guiding-center Hamiltonian dynamics. The guiding- 
center phase-space transformation involves an asymptotic expansion, with a small 
dimensionless parameter 63 = p/Lb <C 1 defined as the ratio of the characteristic 
gyroradius p and the background magnetic- field length scale Lb- This transfor- 
mation is designed to remove the fast gyromotion time scale associated with the 
time-independent background magnetic field Bq associated with an unperturbed 
magnetized plasma 2 . In previous work 9 , this transformation was carried out 
to second order in with the scalar potential $0 ordered at zeroth order in es. 

The results of the guiding-center analysis presented in Ref. 2 are summarized 
as follows. First, the guiding-center transformation yields the following guiding- 
center coordinates (R,p|| , /z, 0, w, t) = Zgc, where R is the guiding-center position, 
P|[ is the guiding-center kinetic momentum parallel to the unperturbed magnetic 
field, /I is the guiding-center magnetic moment, 6 is the gyroangle, and {w, t) are 
the canonically conjugate guiding-center energy-time coordinates (here, time is un- 
affected by the transformation while the guiding-center energy is chosen to be equal 
to the particle energy). Next, the unperturbed guiding-center extended phase-space 
Lagrangian is 

Vgc = -^Al-dR + II {mc/e) d0 - w dt, (40) 

where Ag s Ag -I- (c/e)p||bo is the effective unperturbed vector potential, with 
bg = Bg/i?g; wc, heuccforth, omit displaying the dimensionless guiding-center 
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parameter for simplicity. The unperturbed extended phase-space guiding-center 
Hamihonian is 

Hgc ^ ^ + ^iBo - w = Hg,-w. (41) 

Lastly, from the unperturbed guiding-center phase-space Lagrangian H40() . we ob- 
tain the unperturbed guiding-center Poisson bracket { , given here in terms of 
two arbitrary functions J- and Q on extended guiding-center phase space as [5] 

^ ' ~ mc\d9 dfi dfi de J S*| V dp\i dp\\ 

-^.V^xvg+ f^^-^^V (42) 
eB*|| ^ \dw dt dt dwj' ^ ' 

where Bq = V X Ag and B*^^ = bo • Bq = + (c/e)p||bo • V X bo; note that 
the Jacobian of the guiding-center transformation is J7gc — 'ti-Bqh (i-^-, (Px(Pp = 
Jgc (PRdp^]^ dfi dO). The unperturbed guiding-center Hamiltonian dynamics is, thus, 
expressed in terms of the Hamiltonian H41() and the Poisson bracket (I42|) as = 
{2°', Ti.gc}gc- In particular, the invariance condition /i = for the guiding-center 
magnetic moment follows from the fact that the guiding-center Hamiltonian H41|l is 
independent of the fast gyroangle 6 (to all orders in eg). 

4.2 Perturbed guiding-center Hamiltonian dynamics. We now consider 
how the guiding-center Hamiltonian dynamics is affected by the introduction of 
low- frequency electromagnetic field fluctuations (<i>i,Ai). These fluctuations are 
assumed to satisfy the low- frequency gyrokinetic ordering jlll I15 | : 

e„ = ~ < 1 and p ~ 1, (43) 

where = eBo/mc denotes the charged-particle's gyrofrequency and e^j is a 
small dimensionless ordering parameter associated with the electromagnetic per- 
turbations space-time scales, with lo denoting the characteristic wave frequency, 
fc|| the characteristic parallel wavenumber and kj^ the characteristic perpendicular 
wavevector (both with respect to the unperturbed magnetic fleld Bo). 

Under the electromagnetic perturbations ($i,Ai), the guiding-center phase- 
space Lagrangian 14U|) and Hamiltonian H41|) become 

= Tgco + esTgci and H'g^ = Hgco + esHgci, (44) 

where the zeroth-order guiding-center phase-space Lagrangian Fgco and Hamilton- 
ian Hgco are given by and (|^ . respectively. In what follows, although the 
three small parameters {^Bt ^St ^ui) may be of the same order in practice, we keep 
them separate in order to retain the correct quadratic nonlinearities in the reduced 
Hamiltonian dynamics. In Eq. H44|l . the first-order guiding-center phase-space La- 
grangian Fgci and Hamiltonian Tigci are 

Fgci - ^Ai(R + p).d(R + p) = ^Aig,{-R;fi,e)-d{Ii + p), (45) 

Hgci = e$i(R + p) = e^g,{^R■,^l,^), (46) 

where Aigc(R] fJ-,9) and ^igc(R,; fj.,0) denote perturbation potentials evaluated at 
a particle's position x = R -I- p expressed in terms of the guiding-center position R 
and the gyroangle-dependent gyroradius vector p{fj., 9). Here, the time dependence 
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of ($i,Ai) is not shown explicitly for simplicity and, to lowest order in es, we 
ignore the spatial dependence of p. 

Because of the gyroangle-dependence in the guiding-center perturbation poten- 
tials ($igc, Aigc), the guiding-center magnetic moment ^ is no longer conserved by 
the perturbed guiding-center equations of motion, i.e., fi = 0{es)- To remove the 
gyroangle-dependence from the perturbed guiding-center phase-space Lagrangian 
and Hamiltonian H45l) - (|46|l . we proceed with the time-dependent gyrocenter ph.a,se- 
space transformation 

Z = (R,p||,/i, Z = (R,p||,7Z,^, 

where 2 denote the gyrocenter (gy) extended phase-space coordinates; we note that 
represents a mixed-canonical momentum coordinate (as will be shown later) and 
the time coordinate t is not affected by this transformation. 

The results of the nonlinear Hamiltonian gyrocenter perturbation analysis |15| 
are summarized as follows. To first order in the small amplitude parameter e = es 
and zeroth order in the space-time-scale parameters (ec^, es), this transformation is 
represented in terms of generating vector fields (^1,^2, ■■•) as 

= Z" + eg^ + ••• (47) 

Here, the gyrocenter phase-space Lagrangian is chosen to be of the form 

r = fgc = -AS-dR + {mc/e)jLde - w dt, (48) 

where Aq = Ao + (c/e)pjj bo, so that the gyrocenter Poisson bracket { , j^- = { , }z 
has the same form as the unperturbed guiding-center Poisson bracket H42|l . 

As calculated previously in Eq. (|31|l . the components of the first-order gyro- 
center generating vector field Qi are 

EE Z«} + -^Mgc- {R + P, Z-}, (49) 

where the first-order scalar field Si is determined as follows. The first-order gyro- 
center Hamiltonian is determined from the first-order Lie-transform equation l|34(l 

T-Li = e(^^igc - Aigc--^) - {Si, Tio] = etpigc - {Si, Ho}, 

were V'lgc = ^igc — -^igc • v/c defines an effective first-order perturbation potential. 
The gyroangle-averaged part of this first-order equation yields Tii = e (ipigc), while 
the solution for the scalar field Si is 

g /" — g ^ 

where ipigc = 4'igc — (V'lgc) denotes the gyroangle-dependent part of the first-order 
effective potential V^igc- Next, the second-order term in the gyrocenter Hamiltonian 
is expressed in terms of Eq. as 

where {R -I- p, v} = I/m and {R + p, ^igd = = {R -I- p, Aigd were used. 

Up to second order in the amplitude parameter e, the extended phase-space 
gyrocenter Hamiltonian is, therefore, expressed as 

2/2 2 \ 

n ^Ho + e e(V^ig,) -I- i- (|Ai,,p) - ^ V^i,,}) j , (50) 
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where Tip = p|/2m + fi Bq — w denotes the unperturbed extended guiding-center 
Haniihonian and the gyrocenter paraUel momentum 

6 ^ 

P\\ = P\\ + e-Aigc-bo + 0{ee^,eeB) (51) 

is a mixed-canonical momentum coordinate (i.e., it is kinetic with respect to the 
background vector potential Aq and canonical with respect to the perturbed vector 
potential Ai). Lastly, to first order in e, the low- frequency gyrocenter push- forward 
operator is defined as 

T7'^ = ^ - e (j^{*i5c, ^} + '-A,g,- {R + p, T}^ + 0{e^). (52) 

To lowest order in magnetic- field nonuniformity (i.e., up to e^), the push-forward 
operator is expressed as 

which clearly exhibits the standard nonadiabatic (first term) and adiabatic (second 
and third terms) parts 13 I3U| of the guiding-center Vlasov distribution. 

We have, thus, obtained a reduced (gyroanglc-independent) gyrocenter Hamil- 
tonian description of charged-particle motion in nonuniform magnetized plasmas 
perturbed by low-frequency electromagnetic fiuctuations. At this level, the nonlin- 
ear gyrokinetic Vlasov equation can be used to study the evolution of a distribution 
of test-gyrocenters in the presence of low-frequency electromagnetic fluctuations. 
For a self-consistent treatment that include an electromagnetic field response to 
the gyrocenter Hamiltonian dynamics, a set of low-frequency Maxwell's equations 
with charge and current densities expressed in terms of moments of the gyrocenter 
Vlasov distribution is required. 

5 Variational Principles for Exact and Reduced 
Vlasov-Maxwell Equations 

In this Section, we plan to derive the nonlinear self-consistent gyrokinetic 
Vlasov-Maxwell equations from a reduced variational principle. This reduced vari- 
ational principle will also be used to derive an exact energy conservation law for 
the reduced Vlasov-Maxwell equations. We begin this Section with the variational 
principle for the exact Vlasov-Maxwell equations 31 . 

5.1 Exact Vlasov-Maixwell equations. The variational principle for the 
Vlasov-Maxwell equations is expressed in terms of the action functional [31| 

A ^ - Jd'Z:F{Z)H{Z) + |^F:F, (53) 

where F^i^ = d^A^ — d^Af^ denotes the electromagnetic four-tensor (F = dA is 
an exact two-form); here, we introduced the covariant notation A^ — (— $,A) 
and x^^ = (ci,x), with the Minkowski space-like metric g^^ = diag(— 1, 1, 1, 1). 
We establish the connection between the extended Vlasov distribution T{Z) and 
the time-dependent Vlasov distribution /(z, t) on the six-dimensional phase space 
z = (x, v) by imposing the physical condition w = H{z,t) on the extended phase- 
space Vlasov distribution T{Z\. 

T{Z) = 5[u;-i/(z,i)] /(z,t). (54) 
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Next, the variational principle SA = based on the action functional (|53|l considers 
constrained Eulerian variations for the extended Vlasov distribution J-' defined as 



(55) 



where the Lagrangian variation AT is identically zero for the extended Vlasov 
distribution T and the extended phase-space virtual displacement 8Z^ given by 
Eq. H13(l . is generated by the scalar field S and variations in vector potential bh.: 



SZ" = {Z^ S}^ 



-SA. {x, . 



(56) 



From this definition, the Eulerian variation 
is 



ST = {S, T}z + - (5A.{x, T}z. 



51 of the extended Vlasov distribution 



(57) 



Under Eulerian variations of the electromagnetic potentials (5 A'' = ((5$, &A) and 
the Eulerian variation H57|l of the Vlasov distribution JF, the variation of the action 
functional H53(l . 5 A = J SCd^xdt, can be expressed in terms of the variation of the 
Lagrangian density 



SC = 



d 



— 5A^ F^- 
47r ^ 



d-'v dwS — 



5A, 



1 ^F^''' 



An ^x^^ 



d^vdw {x"", n} 



- / d^vdw S {T, H}z 



T . (58) 



Stationarity of the action functional H53|l with respect to arbitrary virtual phase- 
space displacements generated by S yields the Vlasov equation in extended phase 
space: 

H}z = 0. (59) 
Substituting Eq. into the extended Vlasov equation yields the standard 
Vlasov equation in (6 -I- 1) phase space: 



d dz d , . 
IT. + /(z,i) = 0. 



dt 



dt 9z 



(60) 



d'vv- f{-K,v,t), (61) 



Next, under general variations SA^, of the electromagnetic four-potential, station- 
arity of the action functional H53() yields the Maxwell equations 

— ^ _ 47r- / d^vdw {x", n], T ^ ~ An- 

dxt" cj I ^ iz ^ 

where the w- integration, using Eq. H54(l , yields the standard form of the Maxwell 
equations, with = (c, v). The remaining Maxwell's equations d^^F^^ + d^F^^^ -\- 
daF/^ii, = follow from the fact that F = dA is an exact two-form (i.e., dF — 0). 

5.2 Noether equation for the exact Vlasov-Maixwell equations. Since 
the Vlasov-Maxwell equations (|60(l - (|61|l hold for arbitrary phase-space variations 
generated by S and arbitrary four-potential variations SA^, the variation of the 
Vlasov-Maxwell Lagrangian density H58|l is now expressed as a space-time diver- 
gence known as the Noether equation: 



SC = 



d 
dx^ 



— SA. F^'' 
An ^ 



d'^vdw S — 



H dT 



m dvi, 



(62) 



From the Noether equation, we now derive the energy-momentum conservation law 



T^" = 



(63) 
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for the Vlasov-Maxwell equations (|60() - (|61|l . where T^'' denotes the Vlasov-MaxweU 
energy-momentum tensor. First, we consider arbitrary space-time translation a;^ — > 
a;^ -I- Jx'^ generated by the scalar field S defined as 

S =[mv + {e/c)A]-Sx - w St, (64) 

and potential variations SA^, defined as 

SA^ = f^, Sx"" - d^{A, Sx")- (65) 

Hence, using the spatial components of Eq. Ht)5|) : 

SA = c5<(E + V$) + SxxB - V(A-(5x), (66) 

the space-time translation dx'^ is expressed as 

fe^ = (c6t, (5x) EE {x^", S}z - -SA-{x, x^}z, 

c 



which follows from Eq. I|56|) : note that the phase-space virtual displacements gen- 
erated by Eqs. H64|) and H66|l do not affect the velocity-energy coordinates (i.e., 
<5v = = 6w). Substituting Eqs. and 6C = - d^{6x'' L) into the Noether 

equation H62|) . we find 



F : F 



after the w-integration has been carried out and the terms 



d 



A^Sx" fdFf^" 



+ in^ J (fv v" /(x,v,i) 



47r \ dxf^ 

vanish identically as a result of Maxwell's equations (|61|l . 

The energy conservation law {5x^ = cSt) for the Vlasov-Maxwell equations 
can, thus, be expressed as 

^+V.S-0, (67) 
where the energy density £ and energy-density flux S are 

f(x,t) = i-(|Ep + |Bp) + I d3„||vp/(x,v,i), (68) 

S(x,t) = ^ExB+ /d^i; ^|v|2v/(x,v,i). (69) 
47r J 2 

The momentum conservation law for the Vlasov-Maxwell equations, on the other 
hand, can be expressed as 

dn 

(70) 

where the momentum density 11 and momentum-stress tensor T are 

n(x,i) = + j (fvm^ f{-K,w,t), (71) 

T(x,t) = ^(|Ep |Bp) - ^(EE + BE) 

Svr 47r 

-f J d^v mvv f{x,v,t). (72) 

The exact Vlasov-Maxwell equations (|6U|l - (|61f) have, thus, been derived by a varia- 
tional principle based on the action functional (|53l) involving the extended Vlasov 
distribution and the four-potential A^. Exact energy- momentum conservation 
laws have also been derived by an application of the Noether method. 
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The covariance of the Vlasov part of the action functional (|53ll is used next to 
construct the gyrokinetic Vlasov action functional in which the extended particle 
phase-space coordinates 2, the extended Vlasov distribution JF, and the extended 
particle Hamiltonian 7i are replaced with their gyrocenter equivalents. 



5.3 Nonlinear low-frequency gyrokinetic Vlasov-Maxwell equations. 

The reduced action functional for the low-frequency gyrokinetic Vlasov-Maxwell 
equations [SU |33] is 

^R^ - j d^ZHZ)n{Z) + (|V$|2 - |B|2), (73) 

where we, henceforth, use the notation 

$ EE e(f>i and B = Bq + eVxAi, 

and we omit the overbar to denote gyrocenter coordinates and functions on ex- 
tended gyrocenter phase space. The absence of the inductive part —c~^dtAi of 
the perturbed electric field Ei in the Maxwell part of the reduced action functional 
(I73|l means that the inductive current 9tEi will be absent from Ampere's equation; 
this is consistent with the low-frequency approximation {e^ <C 1) used in nonlinear 
gyrokinetic ordering H43() . 

The variational principle SAr = J SCr d^x = for the nonlinear low-frequency 
gyrokinetic Vlasov-Maxwell equations is based on Eulerian variations for T{Z) 
while variations of the electromagnetic potentials ($, A) are restricted to variations 
of the perturbation potentials <i>i(x, t) and Ai(x, i) only. Variation of Ar with 
respect to 5T{Z) and 5A^'{yi) = ((5$i, 5Ai) yields 



5Ar = - 



J (fZ 5T{Z) n + T{Z) J cPx (^Mi^(x) 



5H 



Mi^(x) 



X 



— ( eV(5$i-V$ - eVX(5Ai-B ). (74) 



Here, the Eulerian variation 5T is constrained to be of the form 

5T = {S, T}, (75) 

where S generates the virtual extended phase-space displacement 6Z'^ = {Z'^, S} 
and { , } is the unperturbed guiding-center Poisson bracket 1421) on extended 
gyrocenter phase space; note that since the gyrocenter momentum coordinates are 
canonical with respect to the magnetic perturbation Ai, the Eulerian variation l|75|l 
has the canonical form. The hmctional derivatives dH/6Aif^{x) in Eq. H74|l . on the 
other hand, are evaluated using the gyrocenter Hamiltonian (|5U|I (to second order 
in e) as 

^-'(?4)V m> 



SAi^ix) 

where Sg^ = (5"^(x — R — p), the push- forward operator T^^ is given to first order 
by Eq. (|5^ . and we used the identity 

/SA'^ 



16 



Alain J. Brizard 



After re-arranging and integrating by parts, the variation H74|) becomes 

+ j (fZ S {T, [w- H)}z + Jd^x{d-J), (77) 

where we have used T{Z) = 6{w — H) /(Z, t) for the gyrocenter Vlasov distribution 
function in extended phase space in the first two terms in SA, while the last term 
in Eq. (|77|l involves the exact space-time divergence 



(78) 

where = {R'^, Ti-Jz denotes the lowest-order gyrocenter four-velocity. Since 
Eq. (|78|l is an exact space-time divergence, it does not contribute to the reduced 
variational principle SAr = 0. 

By requiring that the action functional Ar be stationary with respect to arbi- 
trary variations S and SA'^ (which vanish on the integration boundaries), we find 
the nonlinear gyrokinetic Vlasov equation 

= {T, n}, (79) 

and the gyrokinetic Maxwell equations: the gyrokinetic Poisson equation 

V2$(x) = -47re J Z f (T7I cJ^J = - An e J d' Z {Sl^^j), (80) 

and the gyrokinetic Ampere equation 

VxB(x) = ^ / d'zm {J^'i^SD) EE 1^ jd'Z {.SIJJ). (81) 

If we now substitute J-{Z) = 5{w — H) f{Z,t) into 7{}z — 0, we obtain the 
standard nonlinear gyrokinetic Vlasov equation written explicitly as 

f +f|i5^+ *L.Vi,Vv/-f|i.VH)f ^0. (82, 
ai \Bl^ dp, ei3„|| j \B'^ ) apii 

The nonlinear equations H80|) . (|81|l . and H82I) . with the gyrocenter Hamiltonian (|50|l . 
are the self-consistent nonlinear gyrokinetic Vlasov-Maxwell equations in general 
magnetic field geometry |15| . 

5.4 Gyrokinetic energy conservation law. We now apply the Noether 
method on the gyrokinetic action functional H73|) to derive an exact gyrokinetic 
energy conservation law. By substituting Eqs. O, l|5n|) . and (|ST|l into Eq. (f77|l . 
the variational equation 5A = J 5£ d'^x (we henceforth ignore the subscript R) 
yields the Noether equation 

5C{x) = d-J{x). (83) 

In the Noether method, the variations {S^SA'^^SC) are expressed in terms of gen- 
erators for infinitesimal translations in space or time. 
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Following a translation in time t t + St, the variations S, 6^i, SAi, and SC 
become, respectively, 



s = 


— w5t 


(5$i = 






^5tdtAi 


5C = 


-StdtC 



> . 



(84) 



In Eq. (|84|l . the expression for S satisfies 5t = {i, S}z and the Vlasov-Maxwell 
Lagrangian density is 



C 



1 

8^ 



IBI 



47r 



|V$|' 



EM, 



after the physical constraint 7i = is imposed in the space-time integrand of 
the reduced action functional (|75|l : here, Sem = (|V<i>p + |Bp)/87r denotes the 
electromagnetic-field energy density (in the low- frequency limit). 
By combining Eq. with Eqs. l|75)l and we obtain 
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I 



ii) 



(II) 



{HI) 



dt^i V$ - (E V$ ) X B 
47r An 



(85) 



The term (/) can be written as 



# (- iv^l 



_9 

di 



47r 



y d^ZSH^-R)/ {eT:^<^gc) 



(86) 

where = e^'igc, the identity |V$p = V • ($ V$) - $ V^^ was used, and 
Eq. 1)80(1 was substituted. The term (//) can be written as 



in dt / dt \ 4tt ^ 



V.,_i*.V^|,(87) 



Here, the first term on the right side of Eq. (|87|l cancels the first term in Eq. (|86|l . 
Lastly, the term (///) can be written as 



V - ( V$ X — B 

Att 



= V- 



$ — V X B 

An 



J d'Zf (T-i(v$,,)) 



where Eq. H81|l was substituted. By collecting the remaining terms, we obtain the 
following expression for the local gyrokinetic energy conservation law: 



(89) 
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where the gyrokinetic energy density is 



£{x,t) = J d^ZS^{x-R) f{Z,t) {H ~ 



(90) 



while the gyrokinetic energy density flux is 




(91) 



We obtain the following expression for the global gyrokinetic energy conservation 
law dE/dt — 0, where the global gyrokinetic energy is 



The existence of this exact global energy conservation law allowed the development 
of energy-conserving simulation techniques for nonlinear electrostatic |12| and elec- 
tromagnetic |13l [TSj gyrokinetic equations. 



Through the use of Lie-transform perturbation methods on extended particle 
phase space, we have derived a set of nonlinear low-frequency gyrokinetic Vlasov- 
Maxwell equations describing the reduced Hamiltonian description of gyrocen- 
ter dynamics in a time-independent background magnetic field perturbed by low- 
frequency electromagnetic fluctuations. A self-consistent treatment is obtained 
through a low-frequency gyrokinetic variational principle and an exact gyrokinetic 
energy conservation law is obtained by applying the Noether method. Throughout 
the work, the Lagrangian and Hamiltonian methods (and differential geometry) in 
eight-dimensional extended phase space were emphasized. 

Further developments in gyrokinetic theory not presented here include nonlin- 
ear relativistic gyrokinetic Vlasov-Maxwell equations (34) and nonlinear bounce- 
center Hamiltonian dynamics j35L 136) (based on the asymptotic elimination of the 
fast bounce-motion time scale from the gyrocenter Hamiltonian dynamics). 
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